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We investigate collective spin excitations in two-component fermion condensates with special 
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modes are calculated using Thomas-Fermi theory and the scaling approximation. We demonstrate 
that spin oscillations have more sensitivity to the interaction and the properties of the condensates 
than the density oscillations. 
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I. INTRODUCTION 

Since the realization of the Bose-Einstein condensed (BEG) atomic gases , there has been much 
interest in ultracold trapped atomic systems to study quantum many-body phenomena. Besides the 
Bose-Einstein condensates (BEG) , one can now study degenerate atomic fermi gases ;4] and 

bose-fermi mixtures These systems offer great promise to exhibit new and interesting phenomena 
of quantum many-particle physics. An important diagnostic signal for these systems is the spectrum of 
collective excitations. Such oscillations are common to a variety of many-particle systems and are often 
sensitive to the interaction and the structure of the ground state. 

However, for harmonically trapped gases of either bosonic or fermionic atoms, the frequencies of the 
density oscillations are quite insensitive to the strength of the interaction. In particular, the dipole 
oscillation frequency is completely independent of the interaction As we will see, excitations where 
the two components move out of phase (spin excitations) behave quite differently. 

The underlying theoretical tool to treat dynamic problems of dilute quantum gases is time-dependent 
mean field theory. This reduces to the random-phase approximation (RPA) for small amplitudes, and 
the theory in this form has been applied to density oscillations of these systems [3, 1^ ■ When the single- 
particle spectrum is regular, the long-wavelength excitations are collective and simpler methods can be 
used to calculate the frequencies, in particular with sum rules jO] or the scaling approximation [ijlllj . !!^ . 
In this work we will follow the last approach, which we believe is justified by the extreme regularity of 
the harmonic trapping potential. The scaling method is physically quite transparent, but it is in fact 
equivalent to the theory based on energy- weighted sum rules , as used for example in ref. Q . 

In the dilute limit, the interaction in a two-component fermion condensate is characterized by a single 
number, the scattering length a. We will only consider here the case of positive scattering lengths which 
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correspond to a repulsive interaction. When the interaction is attractive, the system is superfiuid in its 
ground state and the excitation properties are controlled by the energy gap. 

The ground state properties of the two-component fcrmion system with a repulsive interaction was 
treated by Sogo and Yabu [ijj , who showed that there are three regimes, depending on interaction 
strength. For small interaction strengths, the ground state has equal densities of the two components, 
which we call the "paramagnetic" regime. Beyond a certain threshold in interaction strength, the ground 
state becomes "ferromagnetic", that is with a unequal densities of the two components. At very large 
interaction strengths, the minority component vanishes. We shall call this the single-component phase. 
In this work, we will treat systems in the presence of an external magnetic field, which will produce a 
net spin in the ground state even the paramagnetic regime. 

We now briefly describe the spin modes that could be excited by a time-dependent external magnetic 
field. The simplest field to consider is one that is spatially uniform. However, such fields cannot change the 
spatial part of the single-particle wave functions of the condensate and do not induce internal excitations. 
The response to a uniform field is identical to that of a noninteracting ensemble of A'^i — N2 atoms, where 
Ni is the number of atoms in component i. The next field to consider has a dipole spatial dependence. 
Due to the constraints of Maxwell's equations, the overall multipolarity of such a field is J = 2, but 
we shall follow the common terminology calling it dipole. The spin-dipole mode has been discussed 
previously in in refs. 6] and 7|. Ref. 1^ treats the the spin dipole mode as overdamped and does not 
discuss its frequency. In ref. Jj, the spin dipole response was calculated in RPA for a particular value 
of the coupling strength, and it was found to be rather narrow and close to the unperturbed oscillation 
frequency. In our work here, we will calculate the frequency under a variety of conditions: different 
coupling strengths and in the presence of an external static magnetic field. The damping of the mode is 
beyond the scope of this study. 



II. GROUND STATE 



We begin with the expression for the energy of a dilute trapped condensate in the mean-field approxi- 
mation: 

'^''^ [ - 1^ E E + y + ^Irl + nlrlKp, + p,) 

+5P1P2 - 5(m " m) ]■ (1) 

Here ipns are orbital wave functions indexed by spin s = 1,2, ^l,t are the longitudinal and transverse 
frequencies of the trapping field, g is the coupling strength of a contact interaction, /i^ are the magnetic 
moments of the atoms, and B is an external magnetic field. The densities ps are given by the usual sum 
over occupied orbitals: ps ~ 'Yl,°n'^ IV'n.sP- 

We shall use the Thomas-Fermi approximation to evaluate the first term, the kinetic energy. We can 
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then make a change of variables to simpUfy the appearance of the Thomas-Fermi equations, similar to 
the scaling defined in ref. With the scaled variables and the Thomas- Fermi approximation, the 

expression for the energy is, up to an additive constant, 



Ex — 

Here the variables are defined 



A — 1 O 



/m\3 2 
^ = (^) 9^^^' 

^ (3^2)7;j21 T \0) 

The TF equations for the densities rii 2 are derived by variation of the energy (3) with a constraint on 
the total number of particles. This yields 

n\''^ + gn2 = ej - + B (4) 

2/3 , 25 /r-N 

In this equation the Lagrange multiplier e / has the meaning of the Fermi energy. The solution of these 
equations in the absence of a magnetic field is discussed in detail in ref. These authors make an 

additional rescaling to eliminate g, but we do not do that here. For any positive g, there is a value of e/ 
above which the system becomes ferromagnetic. At 5 = 1, the critical point is at e/ = 20/27. Below that 
value, the minimum energy is obtained for the paramagnetic phase, having ni =712. Just above that 
value, both components are present but the densities are unequal. When ey is increased past e/ = 1, the 
system goes into a single-component phase in the center of the trap, at a; = 0. Away from the center of 
the trap, the Fermi energy is effectively reduced, allowing the paramagnetic phase to persist in the outer 
part of the condensate cloud. Thus the minority component will form a hollow sphere, which we will call 
the hollow spin phase. 

Now consider the effect of a magnetic field, li B ^ 0, the densities ni_2 will unequal no matter what 
the value of g. Still, one can distinguish two kinds of behavior near the center of the trap. If the densities 
of both spins decrease as one moves away from the center, we shall call it paramagnetic. However, for a 
certain range of parameters, the density of the minority component may increase initially, moving away 
from the center. This is the ferromagnetic phase. The system can also form the hollow spin phases in the 
presence of a magnetic field. When the external magnetic field B increases more, the minority component 
disappears completely, giving what we call the single-spin phase. 

Density profiles illustrating these three regimes are shown in Fig. ^ The panels show the densities as 
a function of distance a; in a weak magnetic field {g^B = 1.0 x 10^'') and at three different interaction 
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strengths: g = 0.95(a), 1.05(b) and 1.15(c). From top to bottom, the panels show paramagnetic (a), 
ferromagnetic (b), and the hollow spin phases (c). The solid and dashed lines represent the scaled density 
distribution of major and minor components of fermions. 

The phase boundaries as function of B and g are shown in Fig. |5] The upper and lower columns show 
same results, but in the lower column the vertical line is rescaled by factor g^. 

The dashed lines represents the border between the paramagnetic and ferromagnetic phases, which is 
given by d'^ni{x)ldx^ = at a; = 0. 

The solid line represents the results solved by ^2(0) = in cq.©. This line crosses the dashed line 

g — \/2 and B = 1/27. The curvature of the minority component density, d^n2{x) / dx^ at a; = 0, is 
negative when g < 1/V2 and positive when g > 1/V2 This solid line shows the border between the single 
and paramagnetic spin phases when g < 1/^/2, and the border between the ferromagnetic and hollow 
spin phases when g < 1/V2- 

In the hollow phases there are two boundaries for the minority fermion density where n2 = and 
ni =/= 0. On this boundary the density of majority fermion must satisfy 

nf - gni^2B. (6) 

Here we should consider the following function: 

fin) ^ni^gn- 2B. (7) 

The derivative of the above equation is 

^J{n)^nn)^'^-n--^-g. (8) 

Since f{8/27g^) = 0, /(O) = fil/g^) = -2B < and f'{l/g^) < 0, f{n) becomes maximum at 
n = 8/27g^. Thus f{8/27g^) > is required, otherwise f{n) < 0, and there is no boundary of n2, which 
means 77-2 = in all regions. From that we can obtain the condition for B as 

B < B^.. 4 (9^ - 2^) - 2^ " 
The long-dashed represents B = 2/27g^'^ « 0.0745"^ foj. g > 1/^2, which means the border between the 
hollow and single spin phases. 

In Fig. Owe plot the number asymmetry between two components {N1 — N2) /Nt, where Nt = Ni+N2. 
When the external magnetic field is small, the asymmetry of the fermion number is very small and almost 
independent of the coupling in the paramagnetic phase (g < 1.0) except when the coupling is very small, 
where the system is close to the border between the paramagnetic and single spin phases. As the coupling 
increases further and exceeds about g ~ 1.0, the system is changed from the paramagnetic spin phase 
to the ferromagnetic one. The number asymmetry of the fermion number also increases As the external 
magnetic field becomes larger, the change of the number asymmetry is not drastic, though the rapid 
change of the dipole frequency still remains. 
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III. DIPOLE EXCITATIONS 

In this section we derive an expression for the dipole frequency using the scahng method. Without 
loss of generahty, we take the direction of motion along the z-axis. The time-varying density will be 
parameterized with a collective coordinate Ai as [l^ . 

P.{r)= f^P{ri,r2,r^-h), (10) 

where pf\r) is the density distribution of the ground state. Substituting ea. l|lU|) into eq. iQ), we can 
obtain the variation of the total energy up to the order of O(A^) as 

A£;t = Et-eP 

« -]^m^l{\lNi + \lN2) + A{2\^\2-\l-\l), (11) 
where e'^^ is the ground state energy (Ai — 0), and 



(0) 



When we consider the time-dependence of A;, the mass parameter with respect to Ai is obtained as 

Then the classical equation of motion for A is harmonic, giving rise to the following eigenvalue equation 
for the oscillation frequencies, 



m(n| - u'^)Ni -2A 2A W '^i 



■)2 , ,2\ 



2A TO(Oi - uj^)N2 -2A j \\2 

The eigenvalues of this equation are given 



(13) 



c^^^ri? -(^ + ^)±(^ + ^) (14) 
^ ^m,Ni mN2' ^mNi mN2' ^ ' 

One eigenvalue is lu = with the eigenvector Ai = A2. This is the in-phase oscillation of the system 
which follows Kohn's theorem. 
The other eigenvalue is 

a; = + (15) 

rriNi N2 

The eigenvector is given by Xi/Ni — —X2/N2. In this mode, the two components move in opposite 
directions keeping their center of mass at rest. 

Here we should give a comment. If we consider the dipole oscillation in the transverse relation, the 
frequency normalized by the transverse trapped frequency lud/^It is also equivalent to the right-hand-side 
of ea. fT^ . 
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Using the variable transformation explained in the previous sections, the dipole frequency is written as 



In Fig.^we show the frequency of the dipole oscillation in the out-of-phase oscillation as a function of 
the coupling constant, g, with various external magnetic fields, g^B — 1.0 x 10^"' (dashed), 1.0 x 10^^ 
(solid), 1.0x10^^ (long-dashed), and 5.0x10^^ (chain-dotted), respectively. For comparison we also show 
the frequencies in the symmetric system in zero magnetic field (dotted line). 

As the coupling becomes larger, the dipole frequency decreases monotonically until g w 1.0, and sharply 
increases above that. For g < 1.0, where the number asymmetry is small, the frequency lod does not 
strongly depend on the external magnetic field. In the paramagnetic spin phase {g < 1), the density 
distribution of the minor component of fermion n2 is similar to that of the major component ni, and 
the integrand in A > (jlTfl is positive in all regions. As the coupling increases, the size of the fermi gas 
becomes larger, and then the frequency lojj monotonously decreases. 

As the coupling becomes g ^ 1 one component of the fermions is partially converted into the other 
component, and ferromagnetism appears in the central region. In Fig. dn2/dx > for x < Xc 
and dn2/dx < for x > Xc, where Xc ~ 0.32; namely there is a ferromagnetic region for x < Xc- The 
contribution from this ferromagnetic region to A is negative. As the coupling increase, the critical position 
Xc moves to the surface, A becomes smaller, and then ujd increases. 

In the case of a strong magnetic field, the qualitative behavior is similar. In strong coupling g ^ 1 the 
slope of the density function for the minor fermion n2 is smaller in the strong magnetic field than that in 
the weak magnetic field. As the external magnetic field increases, the dipole frequency becomes smaller. 



UJD/n^T) = {l-i(-^ + -^)}5. 



(16) 



where 




(17) 



(18) 
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IV. MONOPOLE OSCILLATION 

In this section we study the monopole oscillation in spherical trap, fir = ^II = ^m- Generally 
the monopole and quadrupole oscillations are coupled, but qualitative properties are not so significant 
difference between the two collective oscillations. In order to study the monopole oscillation, we introduce 
the following scaling: 

Pi{r) = e''^pf{e^'r). (19) 
Under this scaling the total energy becomes 

Et=J2 i^^^'^i + + ^12, (20) 

1=1,2 

where the Ti and Ui arc the kinetic and harmonic oscillator energy parts in the ground state, respectively. 
The interaction energy V12 appears 

^ -gj d\ {Aip(°V^ + A2r^p^°^} + ATuA? + K22\l + 2i^i2AiA2 (21) 

with 

= -f/'i^r{(3 + r|,)pf}{r|,pf}, (22) 

K22 = -f / d'r {r^/^'m + r^^)p'^\ (23) 

^12 = f / d\ {(3 + r|,)p(°)}{(3 + r|:)pf }. (24) 

The mass parameter of the monopole oscillation is given by the mass times the mean square radius. 
Then the monopole oscillation mode and frequencies ojm are obtained by 

[Smw' - Cm] a = 0. (25) 



with 



m 



Xi 



B„ = -| ^J, (26, 





( 


K12 ] 




Cm = 


H 




I K12 


C22 1 





2Ti + 2C/i + ifn ^12 

K12 2T2 + 2U2 + K22 



(27) 



where 



X^ ^ j d^rpiiry. (28) 
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Using the variable transformation explained in the Sec. ITT 



{^?Bm -Cm\\ = (29) 



with 




where 



Bm = \ . , (30) 

\Q X2 ) 

Cm ^ { . . . , (31) 



Xt ^ j (Tx n^{x)x'' (32) 
f. = lj d^x [n,{x)]i, (33) 
i'n^-f/.^.(3...^)(.|^), (34) 



i^22 = ^ {x-^){3n2 + x^), (35) 

= 2ydx(3ni+. — )(3n, + .— )^-yd. [-^^)- (36) 

In the symmetric system pi = the monopole frequencies are given by ajj^j ~ (C\\ + K\2)lmX\ and 
(Cii — Ki2)/mXi. The eigenvector becomes Ai = A2 (in-phase) in the respect to the former frequency, 
and Ai = — A2 (out-of-phase) in the latter frequency. 

In Fig.|Slwe show the frequencies of monopole oscillations as functions of the coupling constant g with 
g^B = 1.0 X 10^"* (a) and 1.0 x 10^^ (b). The solid and dashed lines represent the smaller and larger 
frequencies of two modes, which we call mode-1 and mode-2, respectively. 

In Fig. we also plot the frequencies of out-of-phase and in-phase modes in the symmetric system 
without external magnetic field with the thin long-dashed and dotted lines, respectively. The frequencies 
of the in-phase oscillation in the symmetric system is about ujm ~ 2$!^/, which is the frequency in a 
non-interacting trapped ideal gas, for any asymmetry. As the coupling constant increases, on the other 
hand, the frequency of the out-of-phase mode monotonously decreases. 

In the weak coupling 5 < 1, the frequencies of the mode-1 and mode-2 are almost same as those of 
the in-phase and out-of-phase modes in the symmetric system. As the coupling constant increases, on 
the other hand, the frequency of the mode-2 rapidly increases above g ~ 1. When the coupling constant 
increases further, the frequency of the out-of-phase mode approaches to that of the in-phase mode, and 
the level mixing between the two modes occurs. Even after the level mixing, the frequency of the mode-1 
increases while the frequency of the mode-2 rarely vary. From those we can suppose that the modes- 1 
and -2 are almost out-of-phase and in-phase modes, though the two modes are mixed and exchange their 
roles in strong coupling. 
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In order to confirm that, here, we calculate the mixing angle 9m defined as 

cos6'„i + — ^ smBm, (37) 






where A is an eigenvector obtained from ea. H29|) . When 9„i — 0,Tr and 9m — tt/2, the modes exhibit the 
pure in-phasc and out-of-phase monopole oscillation of two-component fcrmions. 

In Fig.Elwe show the results with the external magnetic fields g^B = l.Ox 10^^ (a) and g^B = l.Ox 10^^ 
(b). In the case of the small magnetic field (Fig. the two modes, the mode-1 and the mode-2, are 
almost clear out-of-phase and in-phase oscillations before the level mixing. Around g ~ 1.2, the modes-1 
and -2 are mixed, and the two collective modes exchange their roles in the strong interaction limit. This 
feature is less clear, but also seen in the case of the strong magnetic field (Fig.lHb)- 



V. SUMMARY 



In this paper we study the spin excitation on the dipole and monopole oscillations in the asymmetric 
two-component fermion condensed system using the scaling method. As for the dipole oscillations the 
in-phase and the out-of-phase modes are completely decoupled even in the asymmetric system. The two 
modes can be coupled in the monopole oscillation, but these two modes are rarely mixed except when 
their frequencies are close. 

In all kinds of oscillations the in-phase motions do not have any special behavior even if the two 
components are largely asymmetric. The frequencies of the in-phase modes are not so different from 
those of the non-interacting system. 

On the other hand the frequencies of these out-of-phase oscillations monotonously decreases as the 
repulsive force between the two kinds of fermions becomes larger. As the coupling further increases, 
the system changes to ferromagnetic and, the frequencies of the out-of-phase oscillations very suddenly 
increase, thought those of the in-phase oscillations do not vary even if the ferromagnetism appears. In 
addition the frequencies and Ml excitation strength show differences between different number asymmetry 
for g > 1, where the ferromagnetism appears. 

Thus significant information of many body systems can be gotten only from the out-of-phase oscilla- 
tions. In this work we describe the two components of fermions with different two spin states. It may 
be thought to be difficult to establish such experiments because we need two kinds of external magnetic 
fields to control the interaction and the asymmetry. However a similar work can be done in systems 
including two different kinds of fermi gases, where we do not need an external magnetic field to construct 
the asymmetry. If we use two kinds of fermions with different masses, we can study a variety of systems. 
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FIG. 1: The scaled density distribution of two components fermion at g = 0.95 (a), g = 1.05 (b) and g = 1.15 (c) 
with with Cf — 20/27 and g^B — 1.0 x 10~*. The solid and dashed lines represents the distribution of the major 
and minor components, respectively. 
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FIG. 2: The phase diagram of the two components of fermion. The symbols "S", "P", "F" and "H" denote 
region of the single-component, paramagnetic, ferromagnetic, single and hollow phases. 
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FIG. 3: The asymmetry of the two components as a function of the coupling constant g., at e/ = 20/27. The 
dashed, sold, long-dashed and dotted lines represent the results with the external magnetic fields ff^B = 1.0x10"*, 
1.0 X 10"^, 1.0 X 10"^ and 5.0 x 10"^, respectively 
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FIG. 4: The frequencies of the dipole oscillations as a function of the coupling constant g. The dashed, sold, 
long-dashed and dotted lines represent the results with the external mag netic fields g^B = 1.0 x 10"*, 1.0 x 10"^ 
1.0 X 10~^ and 5.0 x 10~^, respectively. The dotted line denote the result of the symmetric system. 
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FIG. 5: The frequency of the monopole oscillations, with the external magnetic fields with g'^B = 1.0 x 10"'* 
(a) and 1.0 x 10~^ (b). The thick solid and dashed lines represent the results of the out-of-phase and in-phase 
oscillation modes, respectively. The thin solid and dashed lines denote the result the out-of-phase and in-phase 
oscillation modes of the symmetric system, respectively. 
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FIG. 6: The mixing angle of the two modes of the monopole oscillation with g^B = 1.0 x 10 ■* (a) and 1.0 x 10 ^ 
(b). The thick solid and dashed lines represent the results of the mode-1 and mode-2 , respectively. 



